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Abstract 

We propose a polynomial time /-algorithm (a deterministic algorithm which uses an 
oracle for factoring univariate polynomials over F^) for computing an isomorphism (if 
there is any) of a finite dimensional Fq(x)-algebra A given by structure constants with the 
algebra of n by n matrices with entries from Fq(x). The method is based on computing 
a hnite Fg-subalgebra of A which is the intersection of a maximal Fq[x]-order and a 
maximal i?-order, where R is the subring of Fg(x) consisting of fractions of polynomials 
with denominator having degree not less than that of the numerator. 
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1 Introduction 

Decomposing finite dimensional associative algebras over a field K include the tasks of isolating 
the radical, computing simple components of the radical-free part and finding minimal one¬ 
sided ideals within these simple components. In this paper we consider the case K = ¥q{x) 
where ¥q is the finite field having q elements {q is a prime power). Decomposing algebras over 
Fq(x) can be applied for example to factorization problems in certain skew polynomial rings, see 
the work [9] of Giesbrecht and Zhang and the recent paper [10] by Gomez-Torrecillas, Lobillo 
and Navarro. The first two tasks mentioned above can be accomplished by the polynomial time 
f-algorithm proposed in the work of the hrst and the third authors with Szanto [T7|. The third 
problem. Ending minimal one-sided ideals in simple algebras appears to be more difficult. In 
this paper we propose a solution which works in the special case when the algebra happens to 
be isomorphic to the full matrix algebra Mn{¥q{x)). 
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If ^ is a K-algebra isomorphic to M„(]K), then hnding a minimal left ideal (or, more gen¬ 
erally, finding an irreducible ^-module) is equivalent to constructing an isomorphism (p ■. 
M„(]K). Indeed, the matrices having possibly nonzero entries only in the hrst column form a 
minimal left ideal in M„(]K), so the inverse image under ^ is a minimal left ideal in A. Con¬ 
versely, if M is an irreducible (that is, an n-dimensional) ^-module, then the action of A on 
M gives an isomorphism A = M„(]K). Therefore the task of finding a minimal left ideal is also 
known as the explicit isomorphism problem. 

Recall, that for an algebra A over a field K and a K-basis Oi,... , 0 ^ of A over IK the 
products OiOj can be expressed as linear combinations of the ap. 


QiiOij 'yijlO‘1 T 'yij2(^2 T ■ ■ ■ T 'yijm^^m- 


The elements 'jijk G IK are called structure constants. In this paper an algebra is considered to 
be given as a collection of structure constants. 

Here we consider the explicit isomorphism problem for IK = Fq(x). For the case IK = Fg the 
polynomial time f-algorithm given in ESI by the third author gives a solution. See also [13] for 
related deterministic methods. Recently the first and the third authors with Schicho |T6| found 
an algorithm for solving the explicit isomorphism problem in the case of number fields. Their 
algorithm is a polynomial time ff-algorithm (it is allowed to call oracles for factoring polynomials 
over hnite helds and for factoring integers), assuming that the degree of the matrix algebra and 
the degree of the number field over Q are bounded. They combined algebraic techniques with 
tools from lattice geometry. Some improvements were given in na. Their results have various 
applications, for instance in arithmetic geometry (see 0, 0, 0)- 

The structure of the paper will be the following. First we develop the necessary notions 
concerning polynomial lattices. In Section [2] we summarize the main tools for handling lattices 
over Fg[a;]. The orthogonality defect inequality and the basis reduction algorithm of Lenstra [H] 
are discussed here. We shall also use extensions by Panins [T9] . 

In the next section we state and prove certain facts about maximal orders over polynomial 
rings. Then we use them to construct maximal Fg[a;] and Fg[i] orders in A. The algorithms 
run in polynomial time if one is allowed to call oracles for factoring univariate polynomials over 
hnite helds (it is a polynomial f-algorithm). 

Let R be the subring of Fg(x) consisting of those rational functions where the degree of the 
denominator is at least as high as the degree of the numerator. The main structural result of 
the paper is that the intersection of a maximal i?-order and a maximal Fg[a:]-order is a hnite 
dimensional Fg-algebra which contains a primitive idempotent of A. This theorem and the 
resulting algorithms are described in Section HJ we propose an algorithm to hnd a primitive 
idempotent of A. Finally we arrive at the following theorem: 

Theorem 1. Let A be isomorphic to Mn(¥g{x)), and given by structure constants. Then there 
exists a polynomial (in n and in the size of the structure constants) f-algorithm which finds an 
explicit isomorphism between A and M„(Fq(x)). 

Together with the polynomial time randomized algorithms of Cantor and Zassenhaus [3] 
(or, when g is a power of a prime bounded by a constant, with the deterministic method of 
Berlekamp [2]), this gives a randomized polynomial time solution in general (and a deterministic 
polynomial time algorithm for small characteristic) for the explicit isomorphism problem in the 
special case IK = Fg(a;). We remark that the main ideas of this paper can be extended to the case 
IK = F(a;) for various fields F of constants. (One just needs efficient methods for decomposing 
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finite dimensional algebras over F, and lattice basis redaction over F[a:].) However, extending 
onr algorithms to finding minimal left ideals in algebras which are isomorphic to fnll matrix 
algebras over hnite extensions of Fq(a:) looks more difhcnlt. 

Our main aim was in this work to show the existence of a polynomial time f-algorithm for 
the explicit isomorphism problem over Fq(a;). No attempt has been made to optimize exponents 
and implied constants in the time bound. Those would require substantial further work. Our 
approach, in return, allowed a relatively short description of the methods and arguments. 

2 Lattices over function fields 

Most of our dehnitions and lemmas come from the seminal paper [TB] of A. K. Lenstra. He 
introduced the notion of reduced basis and found an algorithm which hnds a shortest vector 
in polynomial time in lattices over Fq[a:] (he considered sublattices of Fq[a:]™). Note that the 
analogous problem is NP-hard in the case of integer lattices [1]. First we state certain dehnitions 
about Fq[x]-lattices in Fq(x)™'. 

Definition 2. Let f,g & Fq[a:]. Then we set \^\ = deg(/) — deg( 5 f). We will refer to |.| as 
the valuation (or degree^ of an element of¥q{x). We set |0| = —oo. Let v = (ui, ... G 

Fg(x)™. Then the valuation (or degreej of the vector v is |v| = max(|ni|,..., |nm|)- 

Definition 3. L is a full lattice in ¥q{x)”^ if L = {oibi + ■ ■ • + Ombml ca £ where 

bi,..., hm is a basis (over Fq(a ;)) in Fg(a:)™'. 

Definition 4. Let bi, b 2 ,..., b^ G Fg(a:)™. Then the orthogonality defect OD{hi ,..., b^) is 
defined as OD(hi ,..., b^) = 1^*1 “ I det(i?)| where B is the matrix whose columns are 

the bj, (i = 1,... ,m). 

The following lemma is from [18] . However, there it is stated in a slightly weaker form than 
we need it in this paper. So we state and prove the lemma here as well. The proof is also from 

Hi. 

Lemma 5. Let bi,b 2 ,...,bm G Fg(a;)"* be linearly independent and a = where 

ai G Fq[a;]. Then the following holds for every i: 

loil < |a| + OD(bi,... ,bm) - |bi| (1) 

Proof. Consider the a* as unknowns. In this case we have m linear equations and m variables 
so we can use Cramer’s rule. Note that is a basis so the determinant of the coefficient 

matrix B is non-zero. By Cramer’s rule is equal to the quotient of 2 determinants. In other 
words ttj multiplied by the determinant of the lattice is equal to the determinant where the 
Ah column of B is switched to a. Since these two sides are equal, their valuations are equal 
also (on both sides we have elements from Fq(a;)). Note that the valuation of a determinant 
can be bounded from above by the sum of the valuations of its columns. To formalize this last 
sentence: 


IcTjl -|- I det(i?)| < |bi I -|- |b 2 1 |bj_i| -|- |bj_|_i| |bm| -|- |a| 

m 

= |bi| - |bi| + |a|. 

i=l 

After rearranging we obtain the result. □ 
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An implication of this lemma is the following. If we have a vector with small valuation, 
then the coefficients corresponding to a basis are also small, if the orthogonality defect of the 
basis is small. This also suggests that an ideal basis is one whose orthogonality defect is 0. 
This motivates the following definition. 

Definition 6 . A basis bi, b 2 ,..., b^ € Fq(a;)™ is called reduced if OD(hi ,..., b^) = 0. 

Lenstra proposed a polynomial time method [HI Algorithm 1.7] to compute reduced bases 
of sublattices of Fq[x]™'. We quote [HI Proposition 1.14] below. 

Proposition 7. Let bi, b 2 ,..., b^ be over Fq(x) linearly independent vectors from Fq[a:]™' and 
let L be the ¥q[x]-lattice they generate. Let M = maxi<j<m(|bj|) and let M' = max(M, 1). Then 
there exists an algorithm which takes 0{m^M'(0D(hi ,..., b^) + 1)) arithmetic operations in 
Fq and returns a reduced basis Ci,..., of L for which we have |cj| < M (i = 1,... ,m). 

This result can be extended to find a reduced basis of a full lattice in Fq(x)™'. Let us 
assume that we have a basis bi, b 2 ,..., b^ in Fq(x)™'. Let L be the Fg[a:]-lattice generated by 
these vectors and let B be the matrix with columns bi,..., b^- Let 7 be the least common 
multiple of all the denominators of the entries of B. We consider the lattice L' generated by 
7 bi,..., 7 bm. Note that L' G Fq[a:]™. So using Lenstra’s algorithm one can find a reduced 
basis Cl,..., Cm in L'. Note that | det L'\ = \ det L\ +m| 7 |. This implies that choosing b' = ^c, 
we get a reduced basis of L. Since the orthogonality defect of bi,..., b^ is the same as the 
orthogonality defect of ybi,..., yb^, we obtain the following: 

Proposition 8 . Let bi, b 2 ,..., hm be a basis in ¥q(x)"^ and let L be the ¥q[x]-lattice they gener¬ 
ate. Let 7 be the least common multiple of all the denominators for the entries of hi, b 2 ,..., bm- 
Let M = I 7 I + maxi<j<m(|bj|) and let M' = max(M, 1). Then there exists an algorithm which 
takes 0(m^M'(0D(hi ,..., b^) + 1)) arithmetic operations in Fg and returns a reduced basis 
of L. 

Proof. Lenstra’s method m Algorithm 1.7] together with its analysis [T 8 l Proposition 1.14] 
gives the result. □ 

Given an integer k, the set of elements of the lattice whose valuation is smaller than k is 
a finite dimensional Fg-vector space (this is a consequence of Lemma [5j), and a basis of this 
vector space can also be computed efficiently. 

The algorithm of Proposition [ 8 ] finds a reduced basis of a lattice which is given by a basis. 
However, one can ask the following question: what happens if the lattice is only given by an 
Fg[a;]-module generating system? In such situations an algorithm by Panins [HI Algorithm 3.1.] 
is applicable. It finds a reduced basis of a lattice in ¥q(x)^ given by a system of generators. 
We shall make use of the fact that the valuations of the reduced basis obtained by Panins’ 
algorithm will not be greater than those of the given generators. 

3 Maximal orders over Fg[a:] 

3.1 Preliminaries 

In this subsection we assume that i? is a principal ideal domain with quotient field IK and A is 
a central simple algebra isomorphic to M„(]K). Recall that an R-order in ^ is a full R-lattice 
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which is at the same time a subring of A containing the identity element. Maximal orders 
are orders maximal with respect to inclusion. We start with rephrasing |2ni Theorem 21.6], 
specialized to this setting. 

Proposition 9. Let A = HomK(P, V) where V is a vector space of dimension n over K. Let 
L he any full R-lattice in V. Then Homj^(L, L), identified with the subring 

0{L) = {aeA:aL<L} 

of A, is a maximal R-order in A, and all maximal orders are of this form. 

In terms of matrices, the second statement of the theorem gives the following. 

Corollary 10. Assume that A is a maximal R-order in M„(]K). Then there exists an invertible 
matrix P E Mn(K) such that A = PMn{R)P~^ ■ 

Proof. The theorem with with V = K" gives that every maximal R-order in M„(]K) is 0{L) 
for a full i?-lattice L in W^. Let P be a matrix whose columns are an P-basis of L. □ 

We remark that this claim can be found for quaternion algebras in [22l Exercise 4.2]. 

Our eventual aim is to construct a maximal P-order in M„(K). We will construct an initial 
order Aq in a rather straightforward way and iteratively enlarge it. Strictly speaking, our initial 
object Ao will not be an order. We say that an P-subalgebra A of M is an almost R-order in 
A if it is a full P-lattice in A. Thus orders are almost orders containing the identity element 
of A. It turns out that if Aq is an almost P-order, then the P-lattice generated by Aq and the 
identity element of A is an P-order. 

Discriminants enable us to control the depth of chains of (almost) orders and will also be 

useful in representing orders efficiently. The reduced trace, tr a, of an element a of an M is simply 

the trace of a as an n by n matrix. (This is well defined by the Noether-Skolem theorem.) To 

compute reduced traces efficiently, it is not necessary to know an isomorphism A = iip(]K). 

If n is not divisible by the characteristic of K, then tr a is ^ times the trace of the image of 

a under the regular representation of a. In general, the reduced trace can be computed by 

taking an appropriate coefficient of the nth root of the characteristic polynomial of the regular 

representation. This is because the regular representation of A decomposes as a direct sum of 

n copies of the standard n-dimensional (irreducible) representation. 

The bilinear trace form on A is the symmetric bilinear function (a, b) tr ab. As the matrix 

corresponding to an element of an almost P-order A is similar to a matrix with entries of P, 

the reduced trace of any element of A is from P. The discriminant d{A) can be dehned as 

2 

the principal ideal of P generated by the determinant of the Gram matrix (tr 6 j 6 j)T^j^ where 
bi,... ,bn^ are an P-basis for A. It is nonzero and independent of the choice of the basis. We 
can loosely think of d{A) as an element of P, dehned up to a unit of P. As the bilinear trace 
form is non-degenerate, we have the following (see [201 Exercise 10.3]). 

Proposition 11. Let A and T be almost R-orders in A such that ACE. Then (i(r)|(i(A) and 
A = T if and only if d{T) = d{A). 

The following statement gives an P-lattice as an upper bound for P-orders containing a given 
almost order. An extension to more general rings P is used in the proof of [201 Theorem 10.3]. 
For orders over principal ideal domains it is stated explicitly in [151 Proposition 2.2]. As we 
need a slight generalization to almost orders, we give a proof for completeness. 
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Proposition 12. Let A and T be almost R-orders in A such that ACT. Then F C lA where 
d = d{A). 

Proof. Let bi,...,bn^ be an i?-basis for A. Then an element a G F can be written as a = 
_ 2 

with Oj G IK (i = 1,...,n^). For j = 1,... ,n‘^ put fdj = tr abj. Then the elements fdj 
are in R because the elements abj are in the almost order F which is contained in an i?-order 
and hence have reduced trace from R. By linearity, we have '^■aitrbibj = fdj. Cramer’s rule 
gives that each Oj is a quotient of an element of R and d, which means that a G ^A. □ 

An algorithmic consequence is that it is possible to represent i?-orders containing a given 
almost order A as submodules of the factor module |A/A. This will be particularly useful when 
R = Fg[a;], in which case this factor is an deg d-dimensional vector space over ¥q. 

Our algorithm for computing maximal orders is an adaptation of the method proposed by 
the hrst and third authors for the case i? = Z in [T3]. The method is discussed in the context of 
global helds in the Ph. D. thesis of the first author [12]. The algorithm hnds a maximal order 
in a separable algebra over a global held. The algorithm proposed in this paper works also 
for separable algebras in a similar fashion. However, we only consider the case of full matrix 
algebras as some minor details (e.g. those regarding the trace form) are simpler in this case 
(and this is the only case we need later on). 

For completeness, we include proofs of statements that are not rigorously proved for general 
principal ideal rings in [T5] . 

Let M be a full -R-lattice in A. Then the left order of M is dehned by 

Oi{M) = {aeA\aM(Z M}. 

The set Oi{M) is known to be an R-oidei of A, see [201 Chapter 8 ]. It actually follows from 
the fact that Oi{M) is isomorphic to the intersection of two i?-algebras: the image of A under 
the left regular representation and Hom^(M, M) (embedded into HomK(Al, Al)). 

The next two lemmas will be important tools for the algorithm which hnds maximal orders. 
The hrst one reduces the question of enlarging an order over i? to a similar task for /^.^.-orders 
where vr is a prime element of R. Here < K denotes the localization of R at the prime ideal 
Rtt, that is, R^, = {^ : a, (3 & R with vr /[/3}. If F is an i?-order in A, then F^ = i?^F is an 
i?^-order. 

Lemma 13. Let n be a prime element of R and F be an R-order in A. Suppose that J is an 
ideal o/ F^r such that J > ttF^ and Oi{J) > F^. Put I = F n J. Then we have I > tiT and 
Oi{I) > F. 

This lemma is stated for i? = Z in m Lemma 2.7 ]. The proof goes through for any 
principal domain R. We include it for completeness. 


Proof. Clearly / > ttF and / is an ideal of F. We also have J = Rt,I. Let a G Oz(J)\F.n-. Let 
Oi, 02 ,..., at be a generating set of / as an i?-module. Then these elements generate J as an 
i?„-module whence for i = 1 ,..., f we have 


O-iX O-it 

CLQji = ——di 

Pit 


Ai 


( 2 ) 


where aij,{3ij G R and tt does not divide (3ij. Now put (3 = Then f3aai is in I 

{i = 1,... ,t), whence (3al < I and consequently jSa G Oi{I). Finally we observe that (3a is 
not in F since (3 is not divisible by tt and therefore /5a G F would imply a G F^. The proof is 
complete. □ 
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The next simple statement is stated in [151 Proposition 2.8] for R = Z. It enables us to use 
A in place of in computations regarding sufficiently large one or two-sided ideals of A^. 

Proposition 14. Let A be an R-order in A and tt be a prime of R. Then the map 
X -\- G A) induces an isomorphism of rings A/ttA = Aj^/ttA^. 

Proof. Clearly $ : A —)■ A^/ttAtt is an epimorphism of rings. It is straightforward to check that 
its kernel is vrA. □ 

Now we quote some further theorems and dehnitions from na. The next statement is m 
Proposition 3.1]. 

Proposition 15. Let A^ be an RT,-order in A. Then the residue class ring A^ = A.^/7rA^ is an 
algebra with identity element over the residue class field R^, = Rt./'kRt, = R/nR and dimuAl = 
dim;^^A 7 r. //$ : A^ —)■ Att is the canonical epimorphism, then 7rA,r ^ Rad{ATr) = ^~^Rad{AT^) 
and 4) induces a ring isomorphism At,/R ad{AAj = A,,/Rad{A t,) . 

Now we will introduce the important concept of extremal orders: 

Definition 16. Let A^ and Pjr be RT,-orders in A. We say that Ptt radically contains A^ if and 
only ifVT, ^ A^ and Rad{TT,) ^ Rad{AT,). This is a partial ordering on the set of Rt,- orders. 
Orders maximal with respect to this partial ordering are called extremal. 

The next statement is m Proposition 4.1]. 

Proposition 17. An Rt,- order At, of A is extremal if and only if A,, = Oi{Rad{AT,)). 

Finally, we quote m Proposition 4.5]. 

Proposition 18. Let A,, C r,r RT,-orders in A. Suppose that A,, is extremal and P^r is 
minimal among the R^-orders properly containing A^. Then there exists a two-sided ideal I of 
At, minimal among those containing Rad{AT,) such that Oi{I) ft T^r 

3.2 The algorithm 

We start with a high-level description of the algorithm over a general principal ideal domain 
R. Let i? be a principal ideal domain, IK its held of fractions. Suppose that an algebra 
A, isomorphic to M„(]K) is given by structure constants {i,j, k = 1,... ,n^) from IK with 
respect to a basis ai,...,a„2. We assume that these structure constants are represented as 
fractions of pairs of elements from R. Let 5 be a common multiple (e.g., the product or the 
1. c. m.) of the denominators. Then a' = doi {i = 1,... ,n^) will be a basis with structure 
constants G R. Therefore the i?-submodule Aq of A with basis a'l,... ,a '^2 is an almost 
i?-order. 

We shall compute the discriminant d = d{Ao). Let S = {tti, ..., tt,,} be the set of the prime 
factors of d. Observe, that the discriminant of any i?-order conjugate to Mn{R) is 1. This also 
holds for i?,r-orders for any prime element tt. Therefore, by Corollary fTOl and by Proposition fTTl 
Aq^ is a maximal i? 7 r-order for any prime tt not in S. 

Starting with the order A obtained by taking the i?-module generated by Aq and the identity 
element, for each prime in S we test constructively whether A,, is a maximal R,, order using 
the two tests described below. By constructiveness we mean that in the ”no” case we construct 
an i?-order P 3 If any of the tests hnds such a P, then we proceed with P in place of A. 
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Otherwise, if passes the tests for every vr G S' then we conclude that A is already maximal. 
By Proposition ITT] the number of such rounds is at most the number of the prime divisors of d, 
counted with multiplicities. 

The hrst test is used to constructively decide whether A^^ is an extremal i? 7 r-order by checking 
if Oi{Rad{AT,)) = Att fProposition [T7|l . To this end, we compute the ideal I = Rad{AT,) 0 A. By 
Lemma [T3l A passes the test if and only if 0/(1) = A. Otherwise P = 0/(1) is an order strictly 
containing A. To compute I, we work with the n^-dimensional / 2 / 7 ri?-algebra B = A/vrA. From 
Propositions [TT] and [T5] we infer that I is the inverse image of Rad{B) with respect to the 
canonical map A ^ B. 

If Att passes the hrst test, then we proceed with the test of Proposition [181 if there exists 
an ideal J of A.^- minimal among the two-sided ideals properly containing Rad{A„) such that 
Oi{J) > A^, then we construct an i?-order P that properly contains A. Like for the hrst test, 
we can work in the algebra B = A/ttA. Let R,..., Jm denote the minimal two-sided 

ideals of B which contain Rad{B). We have m < iR. Let R denote the inverse image of Jj with 
respect to the map A ^ B. As in the hrst case we obtain, that we have to compute the rings 
Oi{Ii) for z = 1,... ,m. We can stop when A < ORR) is detected, because then we have an 
order properly containing A. 

3.3 The case R = Fg[a:] 

We continue with details of the key ingredients of an efficient algorithm for R = Fq[a:] following 
the lines above. These will give an /-algorithm whose running time is polynomial in the 
size of the input. The input is an array of n® structure constants represented as fractions of 
polynomials. We assume that the numerators are of degree at most djsi and the denominators 
are of degree at most do- Thus the size of the input is around n^{dD + djsf) logg. 

The 1. c. m. of the denominators and hence a basis for the initial order Aq can be computed 
in polynomial time. The degree of this common denominator is at most n^dij, whence Aq will 
have a basis aR , a^ 2 , where each a' is aj, multiplied by a polynomial of degree at most n^do- 
The structure constants for the basis a[,..., aR are polynomials of degree at most n^do + d^- 
The discriminant d = (i(Ao) can be efficiently computed in a direct way following the dehnition. 
The entries of the matrices for the images of a'- at the regular representation, written in terms 
of the basis aR ..., a '^2 are just structure constants for the basis aj,..., a'^ 2 - Therefore these 
entries are polynomials of degree bounded by rA'dD + d^ and hence the entries of the Gram 
matrix of the bilinear trace form are polynomials of degree + 2dAr. To compute d(Ao), let 

n = p^k where p is the characteristic of and k is relatively prime to p. Then the characteristic 
polynomial of a'a', is the nth power of the characteristic polynomial of a'a' as an n by n matrix. 
Therefore it is of the form 

{X^ - (ti + .. .R = {X^^-kRia'iaRX^^-^ + ...y^ 

= X^' -{ktia'.a'jfX^'-P^ + .... 

It follows that d(Ao) is a polynomial Dq of degree at most 2nP'dD + ^inRdN. 

By Proposition [121 we have A < -^Aq for any Fq[x]-order A > Aq. Therefore we can 
represent A as the Fq[x]-submodule A/Aq of the factor module -j|^Ao/Ao. This factor module 
is an n^ deg Do-dimensional vector space over the held F^. In fact, the elements + Aq 
{i = 1 ,..., n^, k = 0,..., deg Do — 1 form an Fg-basis) and we represent A/Ao by an Fg-basis 
written in terms of this basis. Notice that the ideals I whose left order ORI) we compute 


throughout the algorithm are all (left) Ao-submodules of ;^Ao containing -DqAo- Observe next 
that the multiplication of A induces an F^-bilinear map /i from ^Aq/Aq x I/DqI to 
For a G -^A and b E I, one can set 

/i((3 -|- Aq, b -|- DqI^ = Qjb 1. 

This is well dehned as (;^Ao)(-Do-^) = ^ I- Taking an F^-basis bi,... ,bt of I/DqI, the 

factor Oi{I)/ Aq can be computed as the intersection of the kernels of the linear maps fi{-,bi) 
[i = 1,... ,t)- As the dimensions are bounded by polynomials in n and in the degree of Dq, 
for every I possibly occurring in the algorithm, Oi{I) is computable in polynomial time. Given 
an intermediate order A, we can compute the candidate ideals I by computing the radical of 
B = A/gA for the irreducible factors g of Dq and the minimal two-sided ideals of B containing 
the radical and hnally by taking inverse images of these at the map A -E B. As S is an deg g- 
dimensional vector space over F^, its radical and the minimal two-sided ideals containing it can 
be computed in time polynomial in the input size using for example the deterministic method 
of the third author 1211. The minimal two-sided ideals containing the radical, that is, the simple 
components of B/Rad{B) can be found by the deterministic /-algorithm of Friedl and the third 
author [S]. 

For aik e¥g {i = 1,... ,n^, k = 0,... deg Do - 1), the combination XlLi 
of a'l,..., a '^2 has coefficients whose numerators and denominators are polynomials of degree at 
most deg Do < 2n^dD + 2n?dN. Together with a[,..., a'^ 2 , such representatives for an F^-basis 
of A/Aq give a system of generators over Fq[x] for A. When A turns out to be maximal, then 
we can use the lattice reduction algorithm by Panins [IS] to obtain a basis for A consisting of 
combinations of a'^,..., a '^2 with coefficients having numerators and denominators also of degree 
at most deg Dq < 2n^dD + 2n^dN- (Here we make use of the nature of the reduction algorithm: 
it never increases the maximum degree of the coordinates of the intermediate generators.) This 
gives us the following theorem: 

Theorem 19. Let A be isomorphic to Mn{¥g{x)) given by structure constants having numera¬ 
tors and denominators of degree at most dc > 1- A maximal ¥g[x\- order A can be constructed 
by an f-algorithm running in time {n-\- dc + log . The output of the algorithm is an Fg[a;]- 
basis for A whose elements are linear combinations in the original basis of A with coefficients 
which are ratios of polynomials of degree at most (2n® 2n^)dc- 

Notice that Yl^j=Q^A^ — Therefore a fraction of two polynomials in x of 

degree at most d can also be written as a fraction of two polynomials in ^ also of degree at 
most d. Therefore Theorem m gives the following. 

Corollary 20. Let A and dc be as in Theorem, 1X91 Then a maximal ¥g[^]-order A can he 
constructed by an f-algorithm running in time {n-\-dc+ ^ogq)^^^\ The output of the algorithm 
is an Fg[l-]-6aszs for A whose elements are linear combinations in the original basis of A with 
coefficients which are ratios of polynomials (in x) of degree at most (277,® -\- -\- 2n^)dc- 

We remark that we will actually need an Fq[i](-i)-basis for a maximal Fq[i](i)-order. Ob¬ 
viously, for this an Fq[i]-basis for an Fg[i]-order A whose localization at the the prime ^ is 
maximal, will do. Therefore it will be actually sufficient to apply the main steps of the order 
increasing algorithm only for the prime ^ of Fq[i]. 
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4 Finding a rank 1 idempotent in A 

Let R C Fq(x) be the set of rational functions having degree at most 0 (note that the 0 
polynomial has degree —oo hence it also belongs to R). Thus, if /, G Fq[a:], g ^ 0, then ^ ^ R 
iff deg/ < deg^f. It is easy to see that i? is a subring of Fq(a;). Actually R is the valuation 
ring for the valuation — deg of Fq(x). An alternative view is that i? = Fq[i](-i^, the localization 
of the ring Fg[i] at the prime ideal (^). (In fact, one readily verihes that the elements of R 
are precisely the functions of the form /(y)/(?(y), where f,g are univariate polynomials over 
Fq and the constant term of g is not 0.) Thus i? is a discrete valuation ring, and as such, a 
principal ideal ring. 

The main structural result of the paper is the following theorem. It identifies a finite 
subalgebra C of modest size in A, which contains a primitive idempotent of A. 

Theorem 21 . Let A = M„(Fq(a:)) and let K he a maximal Fq[x]-order in A. Also, let R he the 
suhring of¥q{x) discussed above, that is, the set of rational functions of degree at most zero. 
Let A he a maximal R-order in A. Let hi,... , 6^2 he an ¥q[x\-hasis of A, and for j = 1 ,... ,n^ 
let dj he the smallest integer such that -^bj G A. Let dmin = min{dj : 1 < / < dmax = 
ma.x{dj : 1 < j < n^}. Then 

(i) For every element a G A fl A we have a = Yl (^ih, where the ai are polynomials in Fq[x] 
of degree at most n^d^ax - dmin- 

(a) A n A contains a primitive idempotent of A. 

Proof. Let (j) : A ^ Mn(¥q{x)) be an algebra isomorphism such that </(A) = Mn{R). (Such 
a (j) exists by Corollary [TUI) We show that the Fq[a;]-lattice /(A) in M„(Fq(a;)) (the latter 
considered as Fq(a;)"' ) has determinant 1. To see this, let B be the matrix whose columns form 
an Fq[a;]-basis for the Fq[a:] lattice 4>{A)v C Fq(x)” where r is a nonzero vector from Fq(x)”. 
Then 0(A) = BMn{¥q[x])B~^. The claim on the determinant follows from that the standard 
lattice Fq[a:]" has determinant one and from that the conjugation X i-G BXB~^, considered 
as an Fq(a:)-hnear transformation on Fq(a:)"' , has determinant one. For the latter, notice that 
multiplication by B~^ from the right is similar to a block diagonal matrix consisting of n copies 
of B~^, and hence has determinant (deti?“^)"^, while multiplication by B from the left has 
determinant (deti?)”. 

Let C = AnA. AsA = 0“^(M„(i?)), C can be characterized as the set of the elements a 
of A such that 0(a) has no entries of positive degree. As both A and A are F^-algebras, so is 

C. 

Notice that for 0 7 ^ a G .4, the degree of 0(a) G M„(Fq(x)) (the maximum of the degrees 

of the entries of the matrix 0 (a)) is just the minimal (possibly negative) integer r such that 

;40(a) G Mn{R), or, equivalently, x~^a G A. It follows that the degrees of the entries of 4>{bj) 

are bounded by dmax and hence the orthogonality defect of the basis 0(&i),..., 0 (&n 2 ) for 0(A) 

_ 2 _ 

is at most n'^dmax, because |det0(A)| = 0. Therefore, for a = Y^=i^jbj ^ C Lemma [5] gives 
that aj has degree at most n^dmax — dmin, showing statement (i). 

To establish statement (ii), consider an invertible matrix B G Mn{¥q{x)) for which 0(A) = 
B~^Mn{¥q[x])B. Let us consider the lattice Li = B~^¥q[x]'^ in Fq(x)"'. The determinant of 
Li is obviously deti?”^. Let us denote by 6 be the degree of det 5. Let B~^ui,... ,B~^Un, 
with Ui G Fq[a;]"', be an Fq[a:]-basis of orthogonality defect zero for Li. One can obtain such 
a basis by lattice basis reduction. Similarly, let L2 = B'^¥q[x]. Then L2 is an Fq[a;]-lattice 
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having determinant det B. Let B^u[,..., B'^u'^, with u[ G Fq[x]"', be a basis of defect zero for 
L 2 . Now we define a graph. We connect Ui with «'• with an edge if u'^^Ui 7 ^ 0. This dehnes 
a bipartite graph having these 2n vectors as vertices satisfying Hall’s criterion for having a 
perfect matching. (A set of s vectors from ui,... ,Un having less than s neighbors would span 
a subspace of dimension s having an orthocomplement having dimension larger than n — s.) 
By changing the order of UjS we arrange that 7 ^ 0 (i = 1,... ,n). We have 

n n n 

I) = \B-W + E = -^ + ^ = 0, 

j = l j = l j = l 

whence there exists at least one index i, such that the maximum degree of the coordinates of 
B~^Ui and the maximum degree of the coordinates of B^u[ add up to at most zero. Let i be 
such an index and let S resp. S' be the matrix whose first column is Ui resp. m', and whose 
remaining entries are zero. Now Z = B~^SS'^B is a matrix whose entries are of degree at most 
zero. Also, Z G 0(A). Therefore (l)~^{Z) is in C. Furthermore, Z has rank one as it is similar 
to S'S"^ = Uiu'^. Also, as {uiu'^Y = where p = u'-^Ui 7 ^ 0. It follows that the minimal 

polynomial of Z over Fg(a:) as well as that of (l)~^{Z) is —/iX with a nonzero jj, G Fq(a;). As 
(j)~^{Z) G A n A, we have fi & ¥q[x] H R = F^. Now e = is an idempotent in C such 

that 0(e) has rank one. □ 

Remark 22. We give an example of a C* which is not isomorphic to a full matrix algebra over 
Fq(x). Let A = B~^M 2 {¥q[x\)B where B is the following matrix: 



Let T = M 2 {R), i.e. those matrices whose degree is at most 0. Then C = T fl A is generated as 
an ¥q vector space by the following matrices: 

A 0\ fO 0\ fO 0\ fO 0\ fO 0\ 

VO oJ’U oj’^o ij- 

Note that C has dimension 5 over ¥q, hence it cannot be isomorphic to M 2 {¥q). As a matter 
of fact, it is not even semisimple. The radical of C consists of those matrices whose diagonal 
entries are 0. Finally, note that C/Rad C = Fg © Fg. 

For hnding a primitive idempotent of A inside C we can use the method described in the 
proof of the following lemma. 

Lemma 23. Let C he the finite ¥q-algebra from Theoremand let Ci,... ,er be a complete 
system of orthogonal primitive idempotents in C. Then there exists an i such that e* is a rank 
1 idempotent in A. 

Having a basis ofC at hand (a subset of A), one can find such an idempotent by a polynomial 
time f-algorithm. 

Proof. We note hrst, that the identity element of A is in C, hence C has idempotents. Let x G C* 
be an element which is a rank 1 idempotent in A. By Theorem [21] such an x exists. Next observe 
that there exists an index i, for which CiX is not in the radical of C. For, otherwise ^ 

would be in the radical of C, which is impossible, as x is not nilpotent. Let us denote this 
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primitive idempotent e, by e. Since ex is not in the radical of C, the right ideal exC it generates 
in C contains a nonzero idempotent /. Indeed, we can consider this right ideal as an F^-algebra 
which is not nilpotent. Hence if we factor out its radical, then we have a nonzero idempotent 
there ([71 Corollary 2.2.5]), which can be lifted to an idempotent in exC ([71 Corollary 3.1.2]). 
Write / = exy with a suitable y E C. We have ef = e{exy) = e^xy = exy = /. We verify now 
that both fe and e — fe are idempotent elements: 

ifef = /e/e = /(e/)e = ffe = fe 


and 

(e - fef = + {fef - efe - fee = e + fe - fe - fe = e - fe. 

Furthermore, they are orthogonal: 


/e(e - fe) = (fee) - (/e)^ = fe- (/e)^ = 0 


and 

(e - /e)/e = (e/)e - (/e)^ = /e - (/e)^ = 0 . 

Since e is a primitive idempotent, one has either fe = 0 or fe = e. We show that the hrst 
case cannot happen. If /e = 0 then fef = 0. However, fef = P = f which is not zero. This 
implies that fe = e, and e = exye. Since x had rank 1 in e also has rank 1 in A. 

As for the computational part of the statement, hrst one computes a Wedderburn-Malcev 
complement in C: a subalgebra B oi C which is isomorphic to C/Rad{C). This can be done 
in deterministic polynomial time using the algorithm of m Theorem 3.1]. Then we can use 
for example the polynomial time /-algorithms of [ 8 ] and IZB to compute a complete system of 
primitive idempotents in H. To calculate ranks, we can use the fact that for a E A the left ideal 
aA has dimension rn over Fq(x) where r is the rank of a (considered as an n by n matrix). □ 


We prove a bound on and dmax in the case when A and A are the maximal orders 
constructed in Theorem Uni and Corollary [201 respectively. A is an Fq[a:]-order and A is viewed 
as an i?-order here. 


Lemma 24. For the pair of maximal orders as above, we have dmax A (2ra® + 2nf -|- 2'nf)dc 
and drain A —2(2?7,® + + 2n^)dc 

Proof. For short, we write L = {2n^ + + 2n^)dc. Let oi,..., a „2 be the input basis of A 

we use in the algorithms of Theorem [T9| and Corollary [201 We know that the numerators 
and denominators of the structure constants for A are polynomials of degree at most dc. 
Let g*{l/x) be the smallest common denominator of the structure constants when written as 
rational functions in -. The degree of g* is at most n^dc. We know that the g*{l/x)ai are in 
the starting almost Fq[i]-order Aq, hence they are also in A. Also, one can then write 

(-) - P (- 

\x J x^ \x 


where h{y) G Fq[y] and h(0) 7 ^ 0. We have here £ < n^dc. We claim that 
every i. Indeed 


X 


n^df 


-CL{ — 


—£ 

'W 


■9 


X 


CLi. 


e A hold for 

X C 
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Here the first factor is in i?, the second is in A, thus giving the claim. 

We know from Theorem [T2] that every basis element bj of A is a linear combination of the 
Oj with coefficients a* G Fg(x), and the numerator as well as the denominator of a* has degree 
at most L. We claim now that — wr—rbj G A. Indeed, we have 

1 , , _ , 

The hrst factor is in A, the second is in R and the upper bound follows. 

As for drains we observe that the coefficients for the elements of A in the basis {oj} are rational 
functions of degree at least —L (Theorem [12]). Similarly, by Corollary I2II] the coefficients for the 
elements of A in the basis {at} are rational functions of degree at most L. It follows that for 
d < —2L the element -^bj can not be in A, as the coefficient has degree at least L + 1. □ 

Now we turn to the algorithmic task of hnding (an Fg-basis of) C. 

Lemma 25. Let bi,..., bn^ be the ¥q[x]-basis of A constructed by the algorithm of Theorem UR 
and let Ui,... ,m „2 be the R-basis of A constructed by the method of CorollaruWR From these 
data we can construct an Fq-basis of C in deterministic polynomial time. 

Proof. We consider the elements of A as vectors in the basis ui,..., m„ 2 . This way the elements 
of A can be viewed as vectors from Fq(a;)”' in the usual way: an element a & A with a = 
J2'j=i is represented by the vector 

(ai,... ,a„2)'^ G ¥q{xY^. 

Observe, that a vector as above represents an element of A iff |Q;j| <0 holds for every i. 
Consider now the vectors 6' G Fg(x)"' representing the basis elements bi of A. They generate a 
full Fg[a:]-lattice (corresponding to A) in Fg(a;)"' . We next compute a reduced basis ci,..., c „2 
of this lattice. An element 

n2 

a = /SjCj with Pi G Fq[a;] for i = 1,..., 

i=l 

represents an element of C* = A fl A iff |a| <0. We claim that this latter condition is equivalent 
to the set of inequalities 

lAol = lAI + lol < 0, i = l,...,n^. 

Indeed, as the {q} is a reduced Fg[a;]-basis, from Lemma [5] we obtain that 

lAI < |a| + OD{ci ,..., c„ 2 ) - \ci\ = |a| - |cj| (3) 

for every i, hence if |a| < 0 then \PiCi\ < 0 for every i. Conversely, \PiCi\ < 0 for every i 
obviously implies that |a| < 0. We conclude that the elements xRi such that 1 < i < and 
j is a natural number with j + |cj| < 0 form an F^-basis of C. Theorem [21] and Lemma [HI 
provide a polynomial upper bound for the dimension of C over F^, and hence on the number 
of such elements 

The algorithmic subtasks involved here: change of basis from the input basis to the basis 
{ui}, and the lattice basis reduction both can be done in deterministic polynomial time, hence 
from A and A we obtain C in polynomial time. □ 

polynomial bound for the dimension of C follows also simply from the polynomiality of the algorithm 
described here. 
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The main steps of our algorithm for hnding a rank 1 idempotent element e G A are as 
follows. 


1. Construct a maximal Fq[a:]-order A and a maximal i?-order A, by the f-polynomial time 
algorithms of Theorem [19], and Corollary [201 respectively. 

2. Compute an F^-basis of the hnite algebra C = AflA using the polynomial time algorithm 
of Lemma [251 

3. With the polynomial time f-algorithm of Lemma [23] hnd a complete system ei,..., of 
orthogonal primitive idempotents in C, and then select an e* among them which has rank 
1 in A. Finally output this element e = e*. 


Proof of TheoremUi The correctness and the timing for the hrst Step follows immediately from 
Theorem [TO] and Corollary [20] These, and Lemma [M] imply that C admits polynomial size 
description. Then Lemma [251 settles Step 2. Correctness and polynomiality for the last Step is 
provided by Lemma [22] □ 
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